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ABSTRACT 

We have analyzed the low mass X-ray binary (LMXB) candidates in a Chandra observation of the giant 
elliptical galaxy NGC 4472. In a region observed by the Hubble Space Telescope (HST), approximately 
40% of the bright (Lx >10 37 ergs s _1 ) LMXBs are associated with optically identified globular clusters 
(GC). This is significantly higher than the fraction of bright LMXBs in Galactic GCs and confirms that 
GCs are the dominant sites of LMXB formation in early type galaxies. The «4% of NGC 4472 GCs 
hosting bright LMXBs, on the other hand, is remarkably similar to the fraction of GCs with LMXBs in 
every other galaxy. Although statistical tests suggest that the luminosity of a cluster is an important 
driver of LMXB formation in GCs, this appears largely to be a consequence of the greater number of 
stars in bright clusters. The metallicity of GCs is a strong determinant of LMXB specific frequency, 
with metal-rich clusters about 3 times more likely to host LMXBs than metal-poor ones. There are 
weaker dependences on the size of a GC and its distance from the center of the galaxy. The X-ray 
luminosity does not depend significantly on the properties of the host GC. Furthermore, the spatial 
distribution and X-ray luminosity function of LMXBs within and outside GCs are indistinguishable. 
The X-ray luminosity function of both GC-LMXBs and non-GC-LMXBs reveal a break at «3xl0 38 ergs 
s _1 strongly suggesting that the brightest LMXBs are black hole accretors. 

Subject headings: galaxies: general — galaxies:individual(NGC 4472) — galaxies:star clusters — 
globular clusters:general — X-rays:binaries — X-rays:galaxies 



GCs, suggesting that GCs are even more dominant sites 
of LMXB formation in galaxies without recent star forma- 
tion. 

Dynamical processes, such as tidal capture of neutron 
stars in close encounters with other cluster stars (Clark 
1975; Fabian et al. 1975), and interactions between single 
stars and binaries (Hills 1976) in the high density environ- 
ment of GC cores are some of the leading explanations for 
the preponderance of LMXBs in GCs. However, the very 
small number (13) of bright LMXBs in Galactic globulars 
(e.g. Liu, van Paradijs & van den Heuvel 2001; hereafter 
L01]) and other local group GCs has been the biggest im- 
pediment in testing these models in any detail. Thus, in 
order to isolate and understand the primary environmental 
factors driving LMXB formation in GCs it is imperative 
to study them in distant galaxies, especially globular clus- 
ter rich early types. In this Letter we analyze the LMXB 
population of NGC 4472, the giant elliptical in Virgo, and 
its connection to the GC system. 

2. OBSERVATIONS & DATA REDUCTION 

We have analyzed the «40 ks archival Chandra ACIS- 
S3 image of NGC 4472 obtained by Mushotzky et al. on 
June 12, 2000. After standard data pipeline processing, 
as described on the PSU webpage 4 , we identified the X- 
ray point source populations in the 0.5-2.0 keV, 2.0-8.0 
keV, and 0.5-8 keV images. Using WAVDETECT from 
the CIAO 2.2 package, with a threshold of 10~ 6 probability 

1 Based on observations made with the NASA/ESA Hubble Space Telescope, obtained at the Space Telescope Science Institute, which is 
operated by the Association of Universities for Research in Astronomy, Inc., under NASA contract NAS 5-26555, and on observations made 
with the Chandra X-ray Observatory. 

2 Physics & Astronomy Department, Michigan State University, East Lansing, MI 48824. e-mail: akundu, zepf Opa.msu.edu 

3 SISSA/ISAS, via Beirut n. 2-4, 34014, Trieste, Italy, e-mail: maccarone@ap.sissa.it 

4 http:/ /www. astro. psu.edu/xray/acis/recipes/clcan. html 

1 



1. INTRODUCTION 

Early X-ray studies with the Einstein X-ray observa- 
tory revealed that elliptical and SO galaxies are significant 
sources of X-ray emission. While the bulk of the soft X- 
ray luminosity of X-ray bright galaxies can be attributed 
to kT ~ 1 keV thermal emission from hot (10 7 K) gas 
(e.g. Forman, Jones & Tucker 1985), the hard X-ray flux 
is roughly proportional to the optical luminosity, suggest- 
ing a low-mass X-ray binary origin (Trinchieri & Fabbiano 
1985). With the advent of the Chandra X-ray Observa- 
tory, it is now possible to resolve the point source X-ray 
populations likely to be associated with LMXBs in exter- 
nal galaxies such as NGC 4697 (Sarazin, Irwin & Bregman 
2000) and NGC 1399 (Angelini, Loewenstein & Mushotzky 
2001; hereafter A01), and definitively show that LMXBs 
make a significant contribution to the X-ray flux in typical 
ellipticals. Since the stellar populations in such galaxies 
are at least a few Gyrs old (e.g. Trager et al. 2000), con- 
tamination by high mass X-ray binaries is not an issue. 

Globular clusters are especially fertile environments for 
LMXB formation. Even though GCs account for <0.1% 
of the stellar mass in the Galaxy, they harbor about 10% 
of the La->10 36 erg s^ 1 LMXBs (e.g. Verbunt 2002), in- 
dicating a probability of hosting a LMXB that is at least 
two orders of magnitude larger than for field stars. Recent 
Chandra observations indicate that at least 20% of the 
LMXBs in NGC 4697 (Sarazin, Irwin, & Bregman 2001) 
and as much as 70% in NGC 1399 (A01) may reside in 
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of false detections (<1 false source per field) wc detected 
148 total sources. Of these we eliminate four that are 
either within 8" of the nucleus and may be associated ei- 
ther with the central AGN or density enhancements in the 
hot interstellar gas, or have shapes inconsistent with point 
sources. About 10-15 of the 144 candidates are likely to be 
associated with contaminants such as background AGNs 
(Brandt et al. 2000, Mushotzky et al. 2000); the bulk of 
the rest are expected to be LMXBs. The hardness ratios 
(Counts[2-8 keV]) /Counts [0.5-2 keV]) of all the sources are 
consistent with a value of 0.38 at the 2a level and cannot 
be used to distinguish GCs from contaminants, in agree- 
ment with the findings of A01. Wc extract spectra with the 
PSEXTRACT script using all channels from 0.5 - 8 keV. 
Freezing the neutral hydrogen column density of 1.3 x 10 20 
cm~ 2 we fit a power law spectra in XSPEC 11.0 (Arnaud 
1996). A distance of 16 Mpc, a typical value for the Virgo 
cluster (e.g. Macri et al. 1999), is adopted to convert 
the fluxes to luminosity, and for all other calculations in 
this Letter. The faintest LMXB candidate detected has a 
luminosity of 1.1 xlO 37 ergs s" 1 in the 0.5 - 8 keV band 
in which we measure Lx- Detailed reduction procedures 
and source lists will be presented in a companion paper 
(Maccarone, Kundu, & Zepf 2002). 

In addition to the central HST-WFPC2 image of NGC 
4472 previously studied by us (Kundu & Whitmore 2001; 
hereafter KW01) we analyzed 3 other sets of archival 
WFPC2 V & I images of the galaxy and its halo (Fig. 1). 
Using the techniques described in KW01 we identified 825 
GC candidates in the color range 0.8< V-I <1.4 within the 
ACIS-S3 image, excluding the region within 10" of the nu- 
cleus, and measured their luminosity, color and half-light 
radii (r^). Note that we used updated zero points from 
the HST Data Handbook V 3.1 and re-calculated aper- 
ture corrections for all the pointings for this analysis. We 
detect >80% of the total GC population within the ACIS- 
WFPC2 overlap region, with an estimated contamination 
fraction of a few percent (KW01). Given the high den- 
sity of GCs, sub-arcsecond relative astrometry is essential 
to eliminate false matches. We achieve 0.3" r.m.s. rela- 
tive astrometric accuracy by bootstrapping the HST and 
Chandra positions of obvious matches to the ground-based 
Mosaic image of Rhode & Zepf (2001). 

3. RESULTS & DISCUSSION 

Even though the WFPC2 images cover only «20% of 
the ACIS-S3 frame, 72 of the 144 LMXB candidates lie 
within this region. We consider LMXBs to be associated 
with GCs if they are separated by less than 0.65", where 
there is a natural break in the LMXB-GC angular separa- 
tion. Thirty LMXBs satisfy this criterion. We note that 
increasing the matching distance to 1" only adds 1 more 
source, while removing the color criterion adds 2 candi- 
dates with colors inconsistent with GCs. We choose not to 
relax the matching criterion in order to minimise the pos- 
sibility of spurious matches. Thus, 40% of the Ljjf>10 37 
ergs s" 1 LMXBs detected in NGC 4472 lie in GCs. This 
is significantly higher than the 5-15% figure in the Milky 
Way - based on the L01 catalog and Harris (1996) GC dis- 
tances, with the variation due to how transient sources are 
counted - but smaller than the «70% measured in NGC 
1399 (A01). The fraction of LMXBs in GCs and other as- 



sociated numbers reported here represent "snapshot" val- 
ues at the observational epoch. While it is possible that 
there are underlying populations of long-lived transients 
in the field or GCs, the instantaneous values are central 
to addressing scientifically interesting issues such as the 
importance of LMXBs to the X-ray emission from galax- 
ies, their effect on X-ray background measurements and 
the significance of GCs to these values. Hence it appears 
that at present GCs are the dominant sites for LMXB 
formation in early type galaxies, in consonance with the 
suggestion of White, Sarazin & Kulkarni (2002). Thus, it 
is vitally important to study the GC-LMXB connection in 
these galaxies. 

Although the fraction of LMXBs in GCs is different 
across galaxy types, the fraction of GCs hosting LMXBs 
is remarkably similar in all galaxies studied to date. Ap- 
proximately 4% of NGC 4472 GCs host L x >10 37 erg s" 1 
LMXBs, compared to about 4% in NGC 1399 (A01), «2- 
3% in M31 (using DiStefano et al. 2002 and Barmby & 
Huchra 2001), and 1-4% in the Galaxy (using Liu et al. 
2001 and Harris 1996). Thus the formation efficiency of 
LMXBs in GCs must be strongly driven by the properties 
of GCs, rather than that of the host galaxy. 

3.1. Which Clusters Preferentially Form LMXBs? 

The color distribution of the NGC 4472 GCs is clearly 
bimodal (Fig. 2), with a population of blue, metal-poor 
GCs, and one of red, metal-rich GCs, as found in most 
other ellipticals (e.g. Zepf & Ashman 1993, KW01). KMM 
mixture modeling tests (Ashman, Bird & Zepf 1994) re- 
veal a blue peak at V-I=0.98 and a red one at V-I=1.23, 
with a dividing color of V-I=1.10. One of the more strik- 
ing results (Fig. 2) is that there are 3.3 times as many 
LMXBs in red GCs as there are in the blue. Even after 
accounting for the slightly larger number of red GCs, the 
overabundance factor is still «2.7. Previous observations 
of Galactic GCs have also showed such a trend (e.g. Bcl- 
lazzini et al. 1995), although small number statistics make 
it difficult to disentangle it from other GC properties such 
as distance from the center of the Galaxy. 

The top panel of Fig. 2 also suggests a tendency for 
more efficient LMXB formation in GCs in the inner re- 
gion. The lower panel of Fig. 2 clearly shows that LMXBs 
are formed preferentially in the brightest GCs. This effect 
is actually even stronger than indicated when one consid- 
ers that the faint end of the globular cluster luminosity 
function (GCLF) has not been corrected for incomplete- 
ness. Figure 2 also suggests that smaller GCs are favored 
LMXB hosts. Thus, each of these variables could provide 
the physical environment that promotes LMXB formation 
in GCs. In order to better understand the relative im- 
portance of each of these factors we turn to discriminant 
analysis (DA). 

DA is used to weight and combine the discriminating 
variables in such a way that the differences between pre- 
defined groups are maximized (e.g. Antonello & Raffaelli 
1983). Thus each data point is assigned a discriminant 
score of the form F = w\X\ + W2X2 + ...WiXi where F 
is the discriminant score, Wi is the weighting coefficient 
for variable i, and Xi is the i th discriminating variable, 
such that the distribution of discriminant scores of the 
pre-defined groups is maximally separated along the axis 



Kundu et al. 



3 



of this new composite variable. The absolute values of 
the standardized coefficients, Wi, reveal the relative im- 
portance of the associated discriminating variables. In 
certain cases, where the discriminating variables may be 
correlated, the absolute value of the structure coefficients 
- which are the correlations of each variable with the dis- 
criminant function - may give better estimates of the sig- 
nificance of each of the variables. 

Using the SPSSvlO package we performed DA on the 
LMXB and non-LMXB GC populations with V, V-I, r h 
and distance from the center of the galaxy as the vari- 
ables. The standardized coefficients and structure coeffi- 
cients (within brackets) for selected tests are presented in 
Table 1. Since DA is a linear method, we tested the im- 
portance of GC size by using r^, r| and log(r/j) in turn. 
We report only the r/j analysis, which consistently assigns 
the most power to the size parameter. Two random vari- 
ables, a Gaussian, and a uniform distribution are also used 
to gauge the significance of the results. Fig. 3 shows the 
success of DA in separating the two populations. DA of 
the entire sample clearly shows that the luminosity of a 
GC is the most important factor that drives LMXB for- 
mation in GCs. Color and distance appear to be equally 
important secondary factors that drive LMXB formation, 
with size providing a much weaker discriminant. We note 
however that the incompleteness of the faint end of the 
GCLF has a radial dependance that is both a function 
of the background light and differences in the exposure 
times of the WFPC2 images. Since the GCLF is constant 
at all radii (Kundu, Zcpf & Ashman 2002) and the lumi- 
nosity is known to be uncorrelated with color and r^, re- 
stricting the sample to V<23.5 mag GCs, where the com- 
pleteness is «100%, provides a fairer statistical test for 
distance effects. This sample reveals that V, and V-I are 
the two most important discriminating variables, followed 
by smaller contributions from r/, and distance. An inde- 
pendent statistical test using logistic regression analysis 
reveals similar weights for each parameter, strengthening 
our confidence in the results. 

The strong dependence of LMXB specific frequency on 
GC luminosity could simply be a reflection of the larger 
number of stars in more luminous GCs, or it could pos- 
sibly signal more efficient dynamical LMXB formation in 
massive GCs. The inset in Fig. 3 shows that the LMXB 
density per unit GC luminosity is roughly constant with 
V magnitude. Since GC luminosity is proportional to the 
total number of stars for the roughly constant M/L (and 
similar IMFs) in GCs , this suggests that the dependence of 
LMXB frequency on luminosity can largely be attributed 
to the greater number of stars in bright GCs. It is also 
clear from our analysis that metallicity is a significant, 
strong, independent parameter. There are few theoreti- 
cal explanations that can account for this. One possibil- 
ity is the suggestion of Bellazzini et al (1995), that the 
larger radii of metal-rich stars promotes LMXB forma- 
tion in GCs and facilitates mass transfer by Roche lobe 
overflow. While Grindlay (1987) suggested that a similar 
metallicity trend seen in Galactic GCs can be explained 
by a flatter initial mass function in metal-rich GCs, re- 
cent studies suggest that the observed differences in the 
present day mass function depend primarily on GC evo- 
lution, and not metallicity (e.g. Piotto & Zoccali 1999). 



While the rate of LMXB formation depends on the GC en- 
vironment, our observations reveal no obvious trend of Lx 
with GC luminosity, color, size or distance. Further theo- 
retical studies of the effect of GC environment on LMXB 
formation and evolution are required to understand the 
underlying physics. 

3.2. LMXBs in Clusters vs. LMXBs in the Field 

Arc LMXBs in GCs different from those outside GCs? 
Were the non-GC-LMXBs in ellipticals formed in GCs 
and consequently ejected? In order to address some of 
these questions we compare the GC and non-GC popula- 
tions of LMXBs. The mean luminosity of GC-LMXBs, 
< L x >= 2.4 ± 0.6 x 10 38 erg s^ 1 is slightly higher 
than, but statistically indistinguishable from, the non-GC- 
LMXBs at < L x > = 1.8 ± 0.3 x 10 38 erg s" 1 . Simi- 
larly, Kolmogorov-Smirnov tests revealed that the possi- 
bility that the GC and non-GC-LMXBs are drawn from 
the same Lx or radial distance populations cannot be elim- 
inated at greater than the 30% level. However, the exact 
distribution of each of these quantities in the two popu- 
lations may hold vital clues about their origin. A single 
power-law fit of the cumulative luminosity function (Fig. 
4) of each of the LMXB categories can be eliminated at 
greater than 95% confidence level, while a broken power 
law of the form N(> L) = N cut ( T ^) a ^-/>""'^ fits each 
of the populations well. The non-GC-LMXBs and "All" 
source lists were corrected for background/foreground con- 
tamination, estimated from Mushotzky et al (2000), before 
fitting. As is evident from Fig. 4, all populations show 
a break at L CMt «3xl0 38 ergs s _1 . A similar break is ob- 
served in other galaxies; Sarazin et al (2000) argue that the 
"knee" is at the typical Eddington luminosity for spherical 
accretion onto a 1.4 M neutron star, and may indicate the 
break between neutron star and black hole (BH) accretors. 
The luminous GC-LMXB population cannot be explained 
by multiple bright LMXBs in favorable GC environments 
as we see no trend of higher X-ray luminosity in preferred 
GC hosts. Also, the "knee" in the non-GC-LMXB popu- 
lation of NGC 4472 cannot be caused by multiple LMXBs. 
Thus BH-LMXBs appear to be the most plausible expla- 
nation for the break based on present observations. 

If non-GC-LMXBs were formed in clusters they can es- 
cape from GCs by one of two mechanisms, cluster destruc- 
tion - which is more efficient in the inner regions of galaxies 
(e.g. Vespcrini 2000), and should reveal a surfeit of non- 
GC-LMXBs in the inner regions - or dynamical kicks that 
eject LMXBs from GCs. In the latter case, one might ex- 
pect a spatial density distribution at least as extended as 
that of the GCs. If non-GC-LMXBs are formed in the field 
population, one would of course expect the density distri- 
bution to follow the light profile (which is steeper than 
the GC profile). The bottom panel of Fig. 4 suggests that 
the radial density distribution of GC-LMXBs and non-GC- 
LMXBs is roughly similar, and consistent with the GC 
population. At face value this suggests that dynamically 
ejected LMXBs from GCs may account for a significant 
fraction of non-GC-LMXBs. Studies of the spatial prop- 
erties of LMXBs over a larger radial range, especially the 
innermost regions of X-ray faint galaxies, are required to 
probe whether most non-GC-LMXBs in early type galaxies 
have a GC ancestry. Second epoch Chandra observations 
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will further allow discrimination between possible long- 
lived transient LMXBs associated with the field, and the 
more persistent sources that may be expected to form in 
GCs (Piro & Bildsten 2002). 
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Table 1 

Discriminant Analysis Weights 



V V-I Distance Randoml Random2 p-value Correct 

0.89(0.89) -0.41(-0.35) 0.14(0.34) 0.04(0.22) -0.04(-0.07) 0.05(0.10) 0.000 74.3% 

0.89(0.90) -0.41(-0.35) 0.14(0.34) 0.04(0.22) 0.000 74.4% 

BRIGHT GLOBULAR CLUSTERS (V<23.5) 

0.86(0.82) -0.43(-0.40) 0.19(0.23) 0.28(0.30) 0.000 75.5% 



1 Randoml and Random2 are dummy Normal and uniform random variables respectively 

2 The "p-value" denotes the significance of the null hypothesis that there is no discriminating power in the variables based on 
Wilk's lambda statistic, while "Correct" reports the percentage of cases classified in the right group by DA. 
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Fig. 2. — Top: The V-I colors of GCs vs. distance from the center of NGC 4472 and GC color distribution. LMXB-GC matches are 
indicated by filled circles/bins. Bottom: V magnitude of globular clusters vs. half light radius and the GCLF. 
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Fig. 3. — Representative histograms of the discriminant scores of LMXB-GCs and non-LMXB GCs from discriminant analysis using V, V-I, 
distance, and r^ as discriminating variables. The two populations are well separated. Inset shows the LMXB efficiency per unit GC light in 
arbitrary units as a function of V magnitude. 
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Fig. 4. — Top: Histogram of the cumulative X-ray luminosity functions. Note that the histograms have not been corrected for back- 
ground/foreground contaminants. Lines indicate broken power-law fit to the GC-LMXB profile and contamination corrected fits to the 
non-GC-LMXB and all LMXB profiles. A second possible break at fxlO 37 ' 5 ergs s — 1 may be due to incompleteness effects, or reflect the 
counterpart of the break seen in the Galaxy and M31. Bottom: Radial surface density profile of GC and non-GC-LMXBs with Poisson error 
bars. The GC profile is derived from the density of V<23.5 mag GCs (selected so that the completeness is Ril00% everywhere) and shifted 
by an arbitrary amount vertically. The light profile from Rhode & Zepf (2001), shifted by an arbitrary amount, is shown for comparison. 



